# Prep data
xx <- agData_ISAAA_Area %>% filter(Measurement == "Area")
# Plot
mp <- ggplot(xx, aes(x = Year, y = Value / 1000000)) +
  geom_bar(stat = "identity", alpha = 0.25, fill = "darkgreen") +
  geom_line(color = "darkgreen", size = 1.5, alpha = 0.7) +
  scale_x_continuous(breaks = c(1996, 2000, 2005, 2010, 2015),
                     minor_breaks = seq(1996, 2020, by = 1)) +
  theme_agData() +
  labs(title = "Global Area of GE Crops", x = NULL,
       y = "Million Hectares", caption = myCaption)
ggsave("ge_crops_01.png", mp, width = 6, height = 4)

# Prep data
xx <- agData_ISAAA_Area %>% 
  mutate(Value = ifelse(Measurement != "Percent Change", Value / 1000000, Value),
         Unit  = ifelse(Measurement != "Percent Change", paste("Million", Unit), Unit),
         Text  = ifelse((Measurement == "Percent Change" & Year %in% 1997:1999), Value, NA),
         Value = ifelse((Measurement == "Percent Change" & Year %in% 1997:1999), NA, Value))
# Plot
mp <- ggplot(xx, aes(x = Year, y = Value)) + 
  geom_bar(stat = "identity", alpha = 0.5, color = "Black", fill = "darkgreen") + 
  geom_text(aes(label = Text), y = 0, angle = 270, hjust = 1, vjust = 0.5 ) + 
  facet_grid(Measurement + Unit ~ ., scales = "free_y", switch = "y") +
  scale_x_continuous(breaks = c(1996, 2000, 2005, 2010, 2015),
                     minor_breaks = seq(1996, 2020, by = 1)) +
  coord_cartesian(xlim = c(min(xx$Year)+0.4, max(xx$Year)-0.4)) +
  theme_agData(strip.placement = "outside") +
  labs(title = "Global Area of GE Crops", 
       y = NULL, x = NULL, caption = myCaption)
ggsave("ge_crops_02.png", mp, width = 6, height = 4)

# Prep data
xx <- agData_ISAAA_DVDDVG %>% filter(Unit == "Hectares")
xE <- xx %>% top_n(1, Year) %>% pull(Value) / 1000000
# Plot
mp <- ggplot(xx, aes(x = Year, y = Value / 1000000, color = Area)) + 
  geom_line(size = 1.5, alpha = 0.7) + 
  scale_color_manual(name = NULL, values = c("darkgreen", "darkred")) +
  scale_y_continuous(breaks = seq(0, 100, by = 10), 
                     sec.axis = sec_axis(~ ., breaks = xE)) +
  scale_x_continuous(breaks = c(1996, 2000, 2005, 2010, 2015),
                     minor_breaks = seq(1996, 2020, by = 1)) +
  coord_cartesian(xlim = c(min(xx$Year)+0.5, max(xx$Year)-0.9)) +
  theme_agData(legend.position = "bottom") +
  labs(title = "Global Area of GE Crops", x = NULL,
       y = "Million Hectares", caption = myCaption)
ggsave("ge_crops_03.png", mp, width = 6, height = 4)

# Prep data
xx <- agData_ISAAA_DVDDVG %>% filter(Unit == "Percent")
# Plot
mp <- ggplot(xx, aes(x = Year, y = Value, color = Area)) + 
  geom_line(size = 1.5, alpha = 0.7) + 
  scale_color_manual(name = NULL, values = c("darkgreen", "darkred")) +
  scale_x_continuous(breaks = c(1996, 2000, 2005, 2010, 2015),
                     minor_breaks = seq(1996, 2020, by = 1)) +
  theme_agData(legend.position = "bottom") +
  labs(title = "Global Area of GE Crops", x = NULL,
       y = "Percent", caption = myCaption)
ggsave("ge_crops_04.png", mp, width = 6, height = 4)

# Prep Data
xx <- agData_ISAAA_Crop %>% 
  filter(Measurement == "Area", Crop %in% myCrops) %>%
  mutate(Crop = factor(Crop, levels = myCrops))
xE <- xx %>% top_n(1, Year) %>% pull(Value) / 1000000
# Plot
mp <- ggplot(xx, aes(x = Year, y = Value / 1000000, color = Crop)) +
  geom_line(size = 1.5, alpha = 0.7) +
  scale_color_manual(name = NULL, values = myColors_Crops) +
  scale_y_continuous(breaks = seq(0, 100, by = 10), 
                     sec.axis = sec_axis(~ ., breaks = xE)) +
  scale_x_continuous(breaks = c(1996, 2000, 2005, 2010, 2015),
                     minor_breaks = seq(1996, 2020, by = 1)) +
  coord_cartesian(xlim = c(min(xx$Year)+0.5, max(xx$Year)-0.9)) +
  theme_agData(legend.position = "bottom") +
  labs(title = "Global Area of GE Crops", x = NULL,
       y = "Million Hectares", caption = myCaption)
ggsave("ge_crops_05.png", mp, width = 6, height = 4)

# Prep Data
xx <- agData_ISAAA_Crop %>% 
  filter(Measurement == "Percent of GE Area", Crop %in% myCrops) %>%
  mutate(Crop = factor(Crop, levels = myCrops))
xE <- xx %>% top_n(1, Year) %>% pull(Value)
# Plot
mp <- ggplot(xx, aes(x = Year, y = Value, color = Crop)) +
  geom_line(size = 1.5, alpha = 0.7) +
  scale_color_manual(name = NULL, values = myColors_Crops) +
  scale_y_continuous(breaks = seq(0, 60, by = 10), 
                     sec.axis = sec_axis(~ ., breaks = xE, 
                                         labels = paste(xE,"%"))) +
  scale_x_continuous(breaks = c(1996, 2000, 2005, 2010, 2015),
                     minor_breaks = seq(1996, 2020, by = 1)) +
  coord_cartesian(xlim = c(min(xx$Year)+0.5, max(xx$Year)-0.9)) +
  theme_agData(legend.position = "bottom") +
  labs(title = "Percent of GE Area", x = NULL,
       y = "Percent", caption = myCaption)
ggsave("ge_crops_06.png", mp, width = 6, height = 4)

# Prep Data
xx <- agData_ISAAA_Crop %>% 
  filter(Measurement == "Area Change", Crop %in% myCrops) %>%
  mutate(Crop = factor(Crop, levels = myCrops))
# Plot
mp <- ggplot(xx, aes(x = Year, y = Value / 1000000, fill = Crop)) +
  geom_bar(stat = "identity", color = "black", alpha = 0.7) +
  facet_grid(Crop ~ ., scales = "free_y") +
  scale_fill_manual(values = myColors_Crops) +
  scale_x_continuous(breaks = c(1996, 2000, 2005, 2010, 2015),
                     minor_breaks = seq(1996, 2020, by = 1)) +
  coord_cartesian(xlim = c(min(xx$Year)+0.4, max(xx$Year)-0.4)) +
  theme_agData(legend.position = "none") +
  labs(title = "Year to Year Change in GE Crop Hectares", 
       y = "Million Hectares", x = NULL, caption = myCaption)
ggsave("ge_crops_07.png", mp, width = 6, height = 4)

# Prep Data
xx <- agData_ISAAA_Crop %>% 
  filter(Measurement == "Percent Change", Year > 1997, Crop %in% myCrops) %>%
  mutate(Crop = factor(Crop, levels = myCrops))
# Plot
mp <- ggplot(xx, aes(x = Year, y = Value, fill = Crop)) +
  geom_bar(stat = "identity", color = "black", alpha = 0.7) +
  facet_grid(Crop ~ ., scales = "free_y") +
  scale_fill_manual(values = myColors_Crops) +
  scale_x_continuous(breaks = c(1996, 2000, 2005, 2010, 2015),
                     minor_breaks = seq(1996, 2020, by = 1)) +
  coord_cartesian(xlim = c(min(xx$Year)+0.3, max(xx$Year)-0.3)) +
  theme_agData(legend.position = "none") +
  labs(title = "Percent change of GE Crops Year to Year", 
       y = "Percent", x = NULL, caption = myCaption)
ggsave("ge_crops_08.png", mp, width = 6, height = 4)

# Prep Data
xx <- agData_ISAAA_CropPercent %>% 
  filter(Measurement == "Percent", Crop %in% myCrops) %>%
  mutate(Crop = factor(Crop, levels = myCrops))
xE <- xx %>% top_n(1, Year) %>% pull(Value)
# Plot
mp <- ggplot(xx, aes(x = Year, y = Value, color = Crop)) +
  geom_line(size = 1.5, alpha = 0.7) +
  scale_color_manual(name = NULL, values = myColors_Crops) +
  scale_y_continuous(breaks = seq(0, 80, by = 10), 
                     sec.axis = sec_axis(~ ., breaks = xE, 
                                         labels = paste(xE,"%"))) +
  scale_x_continuous(breaks = c(1996, 2000, 2005, 2010, 2015),
                     minor_breaks = seq(1996, 2020, by = 1)) +
  coord_cartesian(xlim = c(min(xx$Year)+0.5, max(xx$Year)-0.7)) +
  theme_agData(legend.position = "bottom") +
  labs(title = "Adoption of GE Crops", 
       subtitle = "Percent of global cropland planted with GE varieties",
       y = "Percent", x = NULL, caption = myCaption)
ggsave("ge_crops_09.png", mp, width = 6, height = 4)

# Prep Data
xx <- agData_ISAAA_CropPercent %>% 
  filter(Crop != "Others", Measurement != "Percent") %>%
  mutate(Measurement = plyr::mapvalues(Measurement, "Area", "Total Area"))
# Plot
mp <- ggplot(xx, aes(x = Year, y = Value / 1000000, color = Measurement)) +
  geom_line(size = 1.5, alpha = 0.7) +
  facet_wrap(Crop ~ ., ncol= 2, scales = "free_y") +
  scale_color_manual(name = NULL, values = c("black", "darkgreen")) +
  scale_x_continuous(breaks = c(2000, 2005, 2010, 2015),
                     minor_breaks = seq(1996, 2020, by = 1)) +
  theme_agData(legend.position = "bottom") +
  labs(title = "Adoption of GE Crops", x = NULL,
       y = "Million Hectares", caption = myCaption)
ggsave("ge_crops_10.png", mp, width = 6, height = 4)

# Prep Data
xx<- agData_ISAAA_Country %>% filter(Measurement == "Area")
yy <- xx %>% filter(Value > 3000000) %>% pull(Area) %>% unique()
xx <- xx %>% filter(Area %in% yy) %>%
  arrange(desc(Value)) %>%
  mutate(Area = factor(Area, levels = unique(.$Area)),
         Value = Value / 1000000)
xE <- xx %>% top_n(1, Year) %>% pull(Value)
# Plot
mp <- ggplot(xx, aes(x = Year, y = Value, color = Area)) + 
  geom_line(size = 1.5, alpha = 0.7) + 
  scale_color_manual(name = NULL, values = myColors_Areas) +
  scale_y_continuous(breaks = seq(0, 80, by = 10), 
                     sec.axis = sec_axis(~ ., breaks = xE[c(1:4,6)])) +
  scale_x_continuous(breaks = c(1996, 2000, 2005, 2010, 2015),
                     minor_breaks = seq(1996, 2020, by = 1)) +
  coord_cartesian(xlim = c(min(xx$Year)+0.5, max(xx$Year)-0.9)) +
  theme_agData() +
  labs(title = "Area of GE Crops", x = NULL,
       y = "Million Hectares", caption = myCaption)
ggsave("ge_crops_11.png", mp, width = 6, height = 4)

# Prep Data
xx<- agData_ISAAA_Country %>%
  filter(Measurement == "Area") %>%
  arrange(desc(Value)) %>%
  mutate(Area = factor(Area, levels = unique(.$Area)))
# Plot
mp <- ggplot(xx, aes(x = Year, y = Value / 1000000)) + 
  geom_bar(stat = "identity", fill = "darkgreen", 
           alpha = 0.6, color = "black") +
  facet_wrap(Area ~ ., scales = "free_y", ncol = 11) +
  theme_agData() +
  theme(axis.text.x = element_text(angle = 90, hjust = 1, vjust = 0.5)) +
  labs(title = "Area of GE Crop by Country", x = NULL,
       y = "Million Hectares", caption = myCaption)
ggsave("ge_crops_12.png", mp, width = 21, height = 8)